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Abstract. 

The stationary nonlinear Schrodinger equation, or Gross-Pitaevskii equation, is 
studied for the cases of a single delta potential and a delta-shell potential. These model 
systems allow analytical solutions, and thus provide useful insight into the features 
of stationary bound, scattering and resonance states of the nonlinear Schrodinger 
equation. For the single delta potential, the influence of the potential strength and 
the nonlinearity is studied as well as the transition from bound to scattering states. 
Furthermore, the properties of resonance states for a repulsive delta-shell potential are 
discussed. 
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1. Introduction 

In the case of low temperatures, the dynamics of a Bose-Einstein condensate can be 
described in a mean-field approach by the nonlinear Schrodinger equation or Gross- 
Pitaevskii equation [1]. We will focus on the one-dimensional case, which can be 
achieved experimentally by a tight confinement in the two other spatial directions (see 
for example [2] and references therein). The nonlinear Schrodinger equation for the 
macroscopic wavefunction is then given by 

+ Vix)+gmx,t)\-^ = m^^, (1) 

where g = Anfi^aN/m is the nonlinear "interaction strength" and N is the number of 
particles in the condensate. The wavefunction is normalized to = 1. In this ansatz, 
one only takes elastic s-wave scattering into account, characterized by the s-wave 
scattering length a. The scattering length a and thus the nonlinearity g are negative for 
an attractive nonlinear interaction and positive for a repulsive one. Another important 
application of the nonlinear Schrodinger equation is the propagation of electromagnetic 
waves in nonlinear media (see, e.g., [3], Ch. 8). 

Analytic solutions of the nonlinear equation for a non- vanishing potential V{x) are 
rare and therefore it is of interest to study such a simple case in some detail. Here 
we study the nonlinear Schrodinger equation for two simple potentials: a single delta 
potential 

V{x) = XS{x), (2) 
modelling a short range interaction, and the delta-shell potential 

{~]~oc for X ^ 

\5ix-a) forx>0 

with a > 0. The delta-shell is a popular model system for the study of resonances 
and decay. We confine ourselves to the stationary case, where the time dependence is 
given by the factor e~''^*/^. Using units with h = 1 and m = 1, the stationary nonlinear 
Schrodinger equation reads 

1 d^ \ 

2 d^ + ^^^^ + 9\Hx)n ^(x) = fi^ix). (4) 

The solutions of equation for the delta potential and the delta-shell potential are 
essentially the ones of the free nonlinear Schrodinger equation. The wavefunction itself 
is continuous, but due to the delta potential, its first derivative is discontinuous at a; = 0, 
resp. X = a: 

lim {^'{a + e) - tp'{a - e)) = 2X^{a). (5) 

One can easily show that this behaviour, well-known for the Schrodinger equation, is 
not changed by the nonlinearity. Furthermore, in the case of the delta-shell potential, 
the boundary condition ip{0) = has to be obeyed. 
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2. Single delta potential 

The single delta potential © is the easiest model for the study of the existence and 
the properties of bound and scattering states. It has been studied in the context of a 
nonlinear flow [4,5], however rather briefly. 

In the linear case, (7 = 0, equation (jlj) with A < supports a single bound state with 
energy Eq = — and a continuous spectrum for E > 0, however, without embedded 
resonances. The normalized bound state wavefunction is 

M^) = V^e^l^l . (6) 
In the following, we will study the modifications of this linear case due to an attractive 
resp. repulsive nonlinearity. By means of the scaling x = x'/\g\, ip = 4''\/\9\^ ^ — ■^'Idl 
and /i = fi'g^ (which conserves the normalization), the parameter (7 in Q can be removed 
up to a sign. Therefore we will fix the nonlinearity to g = ±1 (with the exception of 
section ing . 



2.1. Attractive nonlinearity 

In the case of an attractive nonlinearity, g = —1, the nonlinear Schrodinger equation 
has the well-known bright soliton solution for A = and /i < [6,7], 

iIj{x) = ksech.{k {x — Xq)) with k = a/— 2/i . (7) 

In order to find nonlinear bound states, i.e. normalizable solutions of equation (@]), bright 
soliton solutions of the form (0) for x > and x < are matched at x = by means 
of condition (0). Obviously, the wavefunction iIj{x) has to be symmetric with respect to 
X = and is therefore given by expression ((Tj) for x > and ipi^) = otherwise. 
Inserting this ansatz into equation (0) leads to the condition 

tanh(fcxo) = A/A;. (8) 
Combined with the normalization of the wavefunction, 

\ilj{x)f dx = 2k'^ / sech^(A; (x — xo)) dx 
) Jo 

= 2A;(1 + tanh(/cxo)) , (9) 

this yields 

fc = i-A, i.e. /i = -|(2A-l)2. (10) 
Because of | tanh(A;xo)| < 1, one finds a condition for the existence of a bound state: 

A<Ae=i. (11) 

A bound state exists for any attractive delta potential but also for a repulsive one, 
provided that its strength is not too large. This effect is due to the attractive self 
interaction — |'?/'(x)p which can compensate a limited external repulsion. 

Figure [H shows the wavefunctions for such bound states for three different values 
of the potential strength A. Quite generally, for an attractive delta potential with a 
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Figure 1. Wavefunctions i^ix) of bound states of the nonlinear Schrodinger equation 
(attractive nonlinearity, g = — 1) for a delta potential V{x) = \5{x) with A = —0.2 
( ), A = (— ) and A = +0.2 ( ). 



negative value of Xq, the wavefunction tends to concentrate at the delta potential with 
decreasing A. A repulsive delta potential repels the wavefunction, xq is positive, and one 
observes two peaks of il){x) at a; = ixo that are pushed further away as A is increased 
toward Ac- For A ^ Ac the wavefunction evolves into two infinitely separated bright 
soliton solutions. 

Remarkably, the bound state ceases to exist at a finite negative value of the chemical 
potential 



This difference to the linear equation or the case of repulsive nonlinearity (see below) 
corresponds to the fact that the wavefunction is no longer bound by an external potential 
but by the internal self-interaction. 

For A > Ac, there is no bound state solution any more, but one can actually find 
periodic stationary solutions in terms of Jacobi elliptic functions [6, 7] 



Here L is the period, p G [0, 1] the elliptic modulus of the Jacobi elliptic function cn and 
K{p) denotes the complete elliptic integral of the first kind. The chemical potential is 
related to these parameters by 



These solutions are of course no longer normalizable, and will be denoted as scattering 
states in the following. Such a periodic solution, characterized by three parameters, the 
chemical potential /i, the period L and the shift Xq, has to fulfil only condition (0). Thus, 
for a fixed value of the potential strength A, there exists a variety of solutions, whereby 
the chemical potential p and the period L can be chosen more or less independently. 
The value of xq is then fixed to satisfy condition 




(12) 




(13) 




(14) 
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In the following, we discuss a particular class of solutions that merge continuously 
into the bound state solution when A is decreased below its critical value Ac- Therefore 
we make the ansatz that /i and il){x = 0) depend continuously on the strength A of 
the delta potential at Ac- In fact, we assume the functional relation to be the same for 
A > Ac and A < Ac, i.e. given by equation (fTn|) . and 

'?/'(0) = fcsech (arctanh(A/A;)) , (15) 

respectively. For a given value of A, we construct solutions that fulfil condition 
(0) and yield the desired values of /i and ip{0). Such solutions can indeed be found 
and figure El illustrates such a wavefunction for A = 0.26, just above the critical value 
Ac = 0.25, in comparison to a bound state solution for A = 0.24. 




Figure 2. Wavefunctions of a bound state (— - 
A 0.26) for a repulsive delta potential V{x) 
close to the transition at Ac — 0.25. 



A — 0.24) and a scattering state ( — , 
X6{x) for an attractive nonlinearity 



In the vicinity of the delta potential at a; = 0, both wavefunctions look rather similar 
and thus the transition from a bound to a scattering state seems to be continuous. The 
observed difference between the bound state and the periodic solution for |x| > L 
disappears in the limit A \ Ac because the period L of the Jacobi elliptic solution 
moves toward infinity. 

To explore this transition in some detail, we consider the position xq of the first 
maximum of |'?/'(x)p, as a function of A, given by 



arctanh 



A 



1/2 -A Vl/2-A^ 
for A < Ac and by the solution of the complex equations 



(16) 



^(0) 



k 



P 



A^(0) = k' 



{ kxo{X) 
cn , 

2p-l \y/2p- 1 

( kxo{\) 
sn , 
2p-l \^2p- 1 



p 



and 



/fcxo(A^ 



p 



[17) 



for A > Ac, where k and ipi^) are fixed by equations (fTUI) and (fT3)) as discussed above. 
At Ac, the function Xo(A) shown in figure El on the left has a logarithmic singularity. 
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Figure 3. Transition from a bound to a scattering state. Left: Position xo(A) of the 
first maximum of the wavefunction. Right: Norm per period L. Note that L oo as 
X\Xc- 



On the right of figure El the norm per period j 



L/2 



is displayed, which 

tends to unity at the critical point Ac, i.e. it approaches the bound state normalization. 
Hence the norm is also continuous. 



2.2. Repulsive nonlinearity 

In the case of a repulsive nonlinearity, g = +1, the nonlinear Schrodinger equation has 
the well-known dark soliton solutions for A = [7, 8] : 

ip = y/Jl tax]h{y/Jl{x — Xq)). (18) 

Making such an ansatz separately for x > and x < and matching at x = with 
respect to condition (jS)) yields xq = regardless of the value of A. Remarkably, the 
wavefunction has a zero at x = even for an attractive delta potential. But solutions 
of this kind are of course not normalizable. Another possible solution is 

'?/'(x) = A; cosech(A; (x — Xo)) , with k = a/— 2/i , (19) 

which is usually discarded because of its unphysical singularity at x = Xq. In the case of 
a delta potential, however, this ansatz reveals proper stationary bound states. Assuming 
(jl9|) for X > 0, a short calculation shows that the wavefunction has to be symmetric, 
tpl—x) = iplx). In addition, xq must be negative because otherwise the wavefunction 
would become singular at x = xq. Condition yields 

tanh(fcxo) = k/X (20) 
and the normalization of the wavefunction requires 

r" + 00 

1= / \il!{x)\^ dx = 2k'^ I cosech^(A; (x — Xq)) dx 



|'?/'(x)|^ dx 

oo 

2A;(1 + coth(A;xo)) . 



(21) 



This leads to 
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which must be positive, yielding the condition 

A < Ae = -i , (23) 

i.e. the delta potential must be sufficiently attractive to overcome the repulsive self- 
interaction in order to support a bound state. 




Figure 4. Wavefunctions of bound states of the nonlinear Schrodinger equation 
(repulsive nonlinearity, g — +\) for a delta potential V{x) = X5{x) for three values 

of the potential strength A = — 1 ( ),A = — 0.7(— •— ) and for the critical case 

A = -0.5 (— ) . 



In figure 01 such bound states are displayed for different values of the potential 
strength A. For decreasing values of A, the wavefunction concentrates at the position 
of the delta potential. In the opposite limit, A /" Ac, we see by series expansion of the 
tanh and the sinh functions, that 

Xq ^ Xc = ^-/Xc and — > /ic = (24) 
and that the wavefunction converges to the limiting function 

^c(x) = --^ = (25) 
\x\ — Xc \x\ + 2 

also shown in figure |3] This is in contrast to the case of an attractive nonlinearity where 
the bright soliton peaks move to ±oo at the critical value Ac. 

For A > Ac, one again finds periodic solutions in terms of Jacobi elliptic functions 

|7,8] 

^^,) = ^/E^,MK(pf-^ p), (26) 



L V L 

where L is the periodicity, p G [0, 1] the elliptic modulus of the Jacobi elliptic function 
and K{p) denotes the complete elliptic integral of the first kind. The chemical potential 
is given by 

fx = 8ip+l)KipY/L'. (27) 
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For a fixed value of tlie potential strength A, one again finds a variety of solutions, 
whereby the chemical potential fi and the period L can be chosen more or less 
independently. Note that such periodic solutions can only be found for /i > 0. 

Nevertheless, one can find a lower bound for the period L. From equation ()27p it 
follows that 



For A \ Ac and /i \ the period of the wavefunction L tends to infinity and the 
wavefunction is not periodic any more in this limit. But in this case one cannot find a 
continuous transition to the bound state wavefunction (fT!I|) . For A /" Ac one finds the 
bound state with ip{^) = 1/2 and /i = 0. In contrast, we have i^^O) ^ for \ 
for the periodic solution because of equation (j2Zj). In fact, the elliptic function sn 
evolves continuously into the tanh when the elliptic modulus p tends to unity [9] . 

Actually, there exist Jacobi elliptic functions that merge continuously into the 
cosech as the elliptic modulus p tends to one. These solutions are given in terms of 
the Jacobi elliptic functions ds and cs [9]. But these functions have poles at the zeros 
of the sn and thus are not physical. 



2. 3. Variation of the nonlinearity 

In this section, we will briefly discuss the influence of the mean-field interaction strength, 
i.e. the nonlinearity g, on the solutions of the nonlinear Schrodinger equation for an 
attractive delta potential. We therefore reintroduce the parameter g. 
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Figure 5. Bound state wavefunctions 'ipix) of the nonlinear Sclirodinger equation 
for an attractive delta potential V{x) — —5{x) for different values of the nonlinearity: 
(— ), g = (-), g = +1 ( ) and g = g, = 2 {■ ■ ■) . 

The bound state solutions have already been deduced in the previous sections. In 
figure El the wavefunction of such a bound state is displayed for three different values of 
the nonlinearity g = —1, 0, +1 and a fixed potential strength A = — 1. With increasing 
nonlinearity g, the wavefunction is pushed outward. 
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Figure 6. Dependence of the chemical potential ji ( ) and the parameter xq ( — ) on 

the nonlinearity g < gc = for the bound state of the nonlinear Schrodinger equation 
for an attractive delta potential V{x) = —5{x). 



In both cases of attractive and repulsive nonlinearity, the chemical potential is given 

by 

^ = -l{2\ + g)\ (29) 

which follows directly from the matching condition ^ and the normalization of the 
wavefunction. At a critical value of g, the chemical potential /x becomes zero and the 
bound state ceases to exist. The condition for the existence of a bound state is the same 
as discussed in section 12.21 Reformulated in terms of the nonlinearity parameter g, it 
reads: 

g<g,= -2\. (30) 

When g approaches the critical value gc-, the situation is similar to the case of a 
fixed repulsive nonlinearity g and A Ac as discussed in the previous section. The 
wavefunction at the critical value of g is 

M^) = t4^^- (31) 
\x\ — 2A 

The dependence of xq and /i on the nonlinearity g is illustrated in figure IHl The 
position xo is given by equation (jH} for (7 < and (f20|l for (7 > 0, however with 

k = -\-g/2. 

For g = 0, one finds the well-known value fi = — A^/2, whereas the function ^0(5') 
has a logarithmic singularity. Nevertheless, the bound state wavefunction il){x) evolves 
smoothly into the well-known bound state © of the linear problem for an attractive as 
well as a repulsive nonlinearity. 

For g ^ gc = — 2A, the chemical potential fi tends to zero and Xq tends to the finite 
value 1/A. The disappearance of the bound state if g is increased above gc is similar to 
the effect observed by Moiseyev et al. [10] for a smooth potential V{x) where a bound 
state is transformed into a resonance-like state at a critical nonlinear interaction. 
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3. Delta— shell potential 

In this section we discuss another simple and very popular model system: the delta-shell 
potential. A detailed discussion of the linear three-dimensional delta-shell potential can 
be found in [11]. Here we restrict ourselves to the one-dimensional case 

^^''^^l X6{x-a) forx>0 ^^^^ 

with a > 0. First we briefly resume the basic features of the delta-shell potential in the 
linear case {g = 0), in particular the existence of bound states in an attractive potential 
and resonances in a repulsive one. As we have already discussed the properties of bound 
states in a single delta potential in some detail, we now concentrate on the case of a 
repulsive potential (A > 0). We set h = 1 and m = 1 as above. 

3.1. The linear case 

In the linear case g = 0, the wavefunction in a delta-shell potential is given by 

{sm(kx) for x < a 

■ n ^ 2X . n ^ ■ ^^ . ^^^^ 

sm[kx) + ^ sm(fcaj sm{K[x — a)) for x > a, 

The phase shift 5{k) between incoming and outgoing waves for x > a is easily calculated 
and yields 

^ ,, , cos(2A;a) — 1 . 
The S-matrix S{k) is defined in terms of the the phase shift 5{k) by [12]: 

Bound states correspond to poles of the S-Matrix S{k) on the positive imaginary 
axis. Calculating these poles one arrives at 

e^^'^'^ = 1 - T • (36) 
A 

This equation has a solution on the positive imaginary axis if the condition 

Aa > - ^ (37) 

is fulfilled. This implies that the delta-shell potential has to be sufficiently attractive 
to support a bound state. If the distance a is reduced or A is increased, so that the 
condition (j37jl is not fulfilled any longer, the bound state is lost and one finds a virtual 
state instead. A virtual state corresponds to a pole of the S-matrix S{k) on the negative 
imaginary axis [12]. The wavefunction of such a state diverges exponentially. For 
a — > cxD the delta-shell potential is equivalent to a single delta potential and the energy 
is E ^ -AV2. 
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Naturally there exist no bound states in a repulsive delta-shell potential, but one 
can find resonance states. A resonance is defined by a pole of the S-matrix S{k) in the 
lower half plane [12]. The energy of the n-th resonance is also complex 

£n = kl/2 = En- ir„/2, (38) 

where the imaginary part r„ is interpreted as a decay rate. In the vicinity of a resonance, 
the phase shift 6{k) rapidly changes by an amount of it. 

The amplitude of a resonance wavefunction is enhanced for x < a. This is illustrated 
in figure El for a delta-shell potential of strength A = 10 at a = 1. The ratio of the 
amplitudes on the left-hand side (x < a) and on the right-hand side {x > a) of the 
delta-shell potential, denoted as Ai resp. Ar, is plotted for real values of the energy. 
The peaks of the amplitude ratio Ai/Ar close to the resonances are clearly visible. The 
squared modulus of the wavefunction of the most stable resonance at Si = 4.488 — 0.0631 
is displayed on the right. Nevertheless one has to keep in mind that the wavefunction 
finally diverges exponentially for complex energies S, whereas it is periodic for real 
energies. 




Figure 7. Left: Amplitude ratio Ai/Ar for a delta-shell potential with a = 1 and 
A = 10 in the linear case. Right: Squared modulus of the wavefunction of the most 
stable resonance at £i = 4.488 — 0.063i. 



3.2. Resonances in the nonlinear case 

Now we come back to the nonlinear Schrodinger equation 

1 d2 



^ ^ + X6(x - a) + g\ilj(x)r ] ipix) = fiipix) for x > 0. (39) 
2 dx^ / 

In the following we will only discuss the case of a repulsive delta-shell potential (A > 0). 
By means of a scaling x = x'/s, ip = ip'y/s, A = \'s and /i = /i's^ for s > (which 
conserves the normalization), the number of independent parameters is reduced to two. 
As we are mainly interested in the effects of a varying nonlinearity, the potential is fixed 
by a = 1 and A = 10 in the following examples. In the linear case we find resonances 
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for this potential. Now we want to identify and characterize resonances in the nonhnear 
case as well. 

But the definition of a resonance becomes somewhat ambitious in the nonlinear 
case. A decomposition into incoming and outgoing waves and thus a definition of the 
S-matrix S{k) is not possible. One method widely used to compute resonances in the 
linear case is exterior complex scaling (see e.g. [13]). This technique has also been 
successfully apphed to the nonlinear Schrodinger equation [10,14]. 

We will not adopt this approach but rather look for solutions that can be expressed 
analytically. We have already learned that the real solutions of the free nonlinear 
Schrodinger equation are given in terms of Jacobi elliptic functions. These solutions 
are matched at a; = a to obtain solutions for the delta-shell potential. The chemical 
potential /i of such a solution is real. Thus we can define a resonance neither by a 
complex eigenenergy nor via the S-matrix. In the following we will rather call a state 
a resonance, when its amplitude is resonantly enhanced in the vicinity of the potential, 
i.e. for X < a. 

Let us briefly discuss the time evolution of nonlinear resonances. Note that the 
states 



with a complex chemical potential = /x,. ~ ir/2 do not fulfill the time-dependent 
nonlinear Schrodinger equation, because the norm of these states is not constant. One 
can circumvent this problem by introducing an additional source term or one considers 
the states (j4n|l just as an adiabatic approximation [14]. On the contrary the states 
pnjl with a real chemical potential /i discussed in this paper fulfill the time-dependent 
nonlinear Schrodinger equation but do not decay. 

Furthermore we have to be cautious about the nonlinear parameter g. A meaningful 
definition of the nonlinearity requires that the norm or the amplitude of a solution 
must be fixed in some way, e.g. by \\ip\\ = 1 in section El This is not applicable 
any longer since resonance states are not normalizable. As a global measure of the 
nonlinear interaction we thus define the mean-field potential g\ip{x)\'^, integrated over 
the "interaction-region" < x < a of the external potential: 



3.3. Attractive Nonlinearity 

First we discuss the nonlinear Schrodinger equation with a negative nonlinearity g, 
corresponding to an attractive mean-field interaction. As stated above, the real-valued 
periodic solutions of the free nonlinear Schrodinger equation with a negative nonlinearity 
can be expressed in terms of the Jacobi elliptic function cn [6,7]. Thus, in order to find 
solutions for the delta-shell potential we make an ansatz of the form separately for 




(40) 




(41) 
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X < a and x > a: 



ipiix) = Aicn{AK{pi){f- + \) pi) hrx<a 



^ X + XQ 



Pr I for a; > a. 



The amplitudes Ai^r and the periods Li^j. are given by 

A^rK{pi^r) ^ 8{l-2pi^r)K{p 



\2 



Al,r = ^^^^j^=-^ and f,= ^'^ ~r,^-..l,r) ^ ^^3^ 
\/\g\^l,r ^l,r 

where pi^r are the elhptic parameters of the solution on the left-hand (x < a) and on 
the right-hand (x > a) side of the delta-shell. Clearly one has only one value for the 
chemical potential, whereas the amplitude A, the parameter p and the period L generally 
differ for x < a and x > a. This is different from the linear case, where the period L is 
fixed with the energy. The chemical potential is positive, > 0, if the elliptic parameter 
is restricted to pi^r € [0, 1/2]. 

The boundary condition ip{Q) = is automatically fulfilled by this ansatz. 
Furthermore the wavefunction must be continuous at x = a, whereas its derivative 
is discontinuous according to equation (jSJ, leading to the conditions: 

I. Aicn{ui\pi) =ArCn{Ur\Pr) (44) 

II. 2XAi cn{ui\pi) = - ^^1^^'' sn{ur\pr)dn{ur\pr) + ^^/"^' sn(M;|p;)dn(M;|p;), (45) 

where the abbreviations ui = K{pi){4a/ Li + 1) and Ur = 4:K{pr){a + xq)/ have been 
used. The first condition can be fulfilled by an appropriate choice of xq, as long as 
cn(Mi|pi)| < \Ay.\. Then one can insert the first condition into the second one and 
arrives at 

-cn.{ui\pif -^T^ xo/2 cn{ui\pi)(lii{ui\pi)sii{ui\pi) 



fi l-2pi ' "^'^ ^{l-2piy/^ 

{l-Pr)Pr {1-Pl)pi 



(46) 



(1 - 2p,)2 (1 - 2piy ■ 

As argued above we are looking for solutions whose amplitudes are resonantly 
enhanced for x < a, i.e. for solutions with a maximum amplitude ratio Ai/Ar. This 
ratio is given directly by the elliptic parameters pi^r'- 



A 



Pl{l - 2pr) 



1/2 

(47) 



1 - 2pi)pr 

A resonant enhancement of the amplitude ratio demands that pi ^ pr- 

In order to identify and analyze resonances of the nonlinear Schrodinger equation 
we now calculate the amplitude ratio Ai/A^ as a function of the the chemical potential 
for different values of the effective nonlinearity (^eff- The left-hand side of figure |H1 shows 
the amplitude ratio as a function of yU for an effective nonlinearity g^^i = —5. As in 
the linear case, illustrated in figure resonances can be clearly identified as maxima of 
the amplitude ratio Ai/Aj.. The resonances are, however, shifted to smaller values of fi, 
whereas the width of the resonances remains similar. 
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Figure 8. Amplitude ratio Ai/Ar as a function of the chemical potential ^ for 
an effective nonlinearity ^eff = —5 (left) and for different effective nonlinearities ^eff 
(right). The shift of the resonance marked by an arrow is displayed as a function of 
5off in figure El 



On the right-hand side of figure |H1 the amphtude ratio Ai/Ar{^) is plotted for 
different values of the effective nonlinearity ^^eff- Resonances are clearly identified for all 
values of g^Q, but the shift of the resonance positions is clearly visible in this illustration. 
We note that the resonance heights barely change with g^f^. 

The observed shift of the resonances will be explained in the following. For 
convenience we rather calculate the chemical potential where Ai = Ar at the sides 
of each resonance, in dependence of (/es- These values of the chemical potential will be 
denoted fi^ and fi^ in the following. They are easier to calculate than the resonance 
positions /i„ because pi = pr holds at these values, furthermore this calculation will also 
reveal the influence of g^s on the resonance width. We note that the wavefunction on 
the interval x G [0, 2a] is symmetric (antisymmetric) around x = a for = fi^ (/i = /i^). 

Using both equations PHj) . the chemical potential can be written as 

/i = gA' (l-pi. (48) 



2p 

The elliptic parameter p can be calculated from the relation 

\n\ A'^r'^ 

pK{pf = ^4^. (49) 
Solving this relation for p leads to 

\g\A^L^ 1 flglA^L^y , 



47r2 2 V 47r2 

and inserting this into equation ()48|1 . we flnd the desired dependence of the chemical 
potential on the nonlinear interaction 

= 1 + + 0^9^)) ■ (51) 



L2 V 87^2 / 

Formula is valid for both /x^ and /x^. Now we insert the speciflc values of the 
period L and replace gA^ by the effective nonlinearity ^^efr- At /z^ the period of the 
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wavefunction is = 2a/n, i.e. ipi^^) = 0- Equation (j4ip for the effective nonlinearity 
can be easily evaluated in lowest order in p, since then the elliptic function cn equals a 
cosine, which yields Qcs ~ gA^a/2. This finally leads to 



'2^ 



+ 



2a 



(52) 



Similarly one obtains an expression for p,^. In the linear case the period is given 



by the solution of the implicit equation 
'27ra\ 2iT 



tan 



r < 



XL< 



(53) 



For the example illustrated in figure El (a = 1, A = 10 and n = 3) one has = 0.7215 . 
The change of with Qcs is negligible. Again equation ()41|) for the effective nonlinearity 
is readily evaluated in lowest order in p and yields 

gaA^ f ^ sin(47ra/L<)^ gaA^ ( ^ 1 



fl'eflf 



47ra/L< 



Aa(l + (2vr)V(AL<)2) 



Inserting into equation (j5ip . one finally arrives at 



27r2 



+ 



3fi'efr 



sin(47ra/L<; 



(54) 



(55) 



(L<)2 ' 2a y A7:alL< 
The same results are obtained in the case of a repulsive interaction ((7 > 0, see 
below). Thus we compare the approximations ()52|1 and ()55j) to the numerically exact 
results for g <Q and g > Q together in figurelHl We considered the resonance with n = 3, 
that is marked with an arrow in the figures |H1 and El We observe a good agreement. 
Furthermore the positions /in=3 of the resonances are displayed. 
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Figure 9. Left: Dependence of the resonance position /i„ (o) and /ij-J' (+) and /i^ 
(x) on the effective nonhnearity ^off for n — i. The sohd hnes are the approximations 
(|5^ and lf55|l . Right: Width of the resonance, defined as A/i = /Xg — /i^ (+) and as 
FWHM (o). 

From the different scaling of p^ and p'^ we conclude that the resonance width 
also changes with the effective nonlinearity. In fact, the width increases almost linearly 
with (yfeff and the resonances become slightly asymmetric. The dependence of the width 
A/i = yUg — yUg on the effective nonlinearity g^s is illustrated in figure El on the right. 
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It should not be concealed that also bound stated can exist in a repulsive delta- 
shell potential due to the attractive self-interaction, falling of as sech(A/— 2/i(x — xq)) 
for X > a. However we will not consider these states here as we already discussed a 
similar phenomenon for the single delta potential. 



3.4- Repulsive Nonlinearity 



As stated above, the real non-singular solutions of the free nonlinear Schrodinger 
equation with a repulsive nonlinearity can be expressed in terms of the Jacobi elliptic 
function sn [7,8]. Thus we make the ansatz: 



ipi^x) 




AisnUKipi)f- 



Pi 



ArSn{AK{pr)^^ 



P1 



for X < a 



for X > a. 



The amplitudes Ai^r and the periods Li^,. are now given by 



A 



A^/p{^K{pi^,) 



Lr 



\9\L 



and fi 



8{pi,r + l)K{p 



Lr) 



Lr 



^Lr 



(56) 



(57) 



where pi^r ^ [0, 1] are the elliptic parameters of the solution on the left {x < a) and on 
the right {x > a) of the delta-shell. 

The boundary condition tp{0) =0 is automatically fulfilled by the ansatz (j^Uj) . The 
remaining conditions for the wavefunction and its derivative at a; = a (cf. equation (jS))) 
read: 



I. Aisn{ui\pi) 

II. 2XAisn{ui\pi 



Ar Sn{Ur\Pr) 

^^''^''-cn{ur\Pr)dn{ur\Pr) " ^^^cn{ui\pi)dn{ui\pi 



(58) 
(59) 



Lr Li 

where the abbreviations ui = 4:K{pi)a/Li and Ur = 4:K{pr){a + xo)/Lr have been used. 

If |A/sn(M;|p;)| < \Ar\ the first condition can always be fulfilled by an appropriate 
choice of the "phase shift" Xq. Inserting the first condition into the second one and using 
the addition theorems of the Jacobi elliptic functions one arrives at 



/i + 1 

Pr 



sn {ui\pi) + 



{pi + 1)3/2 ^ 



cn{ui\pi)dn{ui\pi)sn{ui\pi 



Pi 



(60) 



(61) 



{pr + 1)2 {pi + 1)2 

The amplitude ratio Ai/Ar is given by 

Pl{Pr + 1) 1 ^''^ 
_{pi + l)pr_ 

in terms of the elliptic parameters. A resonant enhancement of the amplitude ratio 
demands that pi ^ pr- 

Again we calculated the amplitude ratio Ai/Ar as a function of the chemical 
potential /x for different values of the effective nonlinearity (^efr- The results are illustrated 
in figure El The left-hand side shows the amplitude ratio Ai/A^ for an effective 
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Figure 10. Left: Amplitude ratio as a function of the chemical potential /i for an 
effective nonlinearity (/off = 5- The condition H()3|l for the amplitude ratio is plotted 
as a dashed line. Right: Amplitude ratio as a function of the chemical potential fjL for 
different effective nonlinearities. The shift of the resonance marked with an arrow is 
displayed as a function of ^cff in figure El 



nonlinearity ^^eff = 5, what should be compared to figure [3 and figure |H1 The first 
observation is that one cannot find solutions for all values of /i. In fact there exist no 
solutions with an amplitude ratio below a certain threshold. Resonances are still clearly 
identified as maxima of the amplitude ratio. Again the resonance positions are shifted 
in comparison to the linear case. 

On the right-hand side the amplitude ratio is plotted for different values of g^a- 
One observes that the solutions cease to exist with an increasing effective nonlinearity, 
whereas the resonances survive longest. The resonances are shifted similarly to the case 
of an attractive interaction. 

Solutions with a small amplitude ratio Ai/A^ cease to exist when ^feff is increased. 
In fact, condition cannot be fulfilled any longer if the amplitude ratio Ai/A^. drops 
below a certain threshold. A condition for the existence of a solution can be derived 
from the equations (fHTjl and (jKTj) and yields 

iiy> ji^.M. (62) 

ArJ Pl + l H 

Inserting ^fes ~ gC'A^/2 on the right hand side, one is led to the approximation 

4l>fi*£V'\ (63) 

As a consequence solutions apart of the resonances with small amplitude ratios cease 
to exist when g^w is increased. This approximate condition is well confirmed by the 
numerical exact results displayed in figure Cni 

The shift of the resonances is understood in the same way as in the case of an 
attractive interaction. The chemical potential is now given by 

+ (64) 
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while equation (j^UI) still holds. Inserting into equation and expanding up to the 
linear term in gA^ again leads to equation ()51|) . Thus one arrives at the same results 
as in the case of an repulsive interaction, in particular at the equation (jHSj) for /x^ and 
equation (j^S)) for fi^. The results forn = 3 are displayed in figurelHl The approximations 
agree well with the numerical exact results. From the different scaling of /i^ and fi^ 
with g^s we conclude that a repulsive nonlinearity increases the resonance width. 

4. Conclusion 

In this paper we analyzed the properties of bound, scattering and resonance states of 
the nonlinear Schrodinger equation using two simple model systems. 

Bound, i.e. normalizable, states were calculated and analyzed for a single delta 
potential. New features occur in the case of an attractive nonlinearity, as states are 
no longer bound by an external potential but by the internal interaction. In this case 
bound states can exist despite a repulsive external potential, and they cease to exist at 
a negative value of the chemical potential. In addition we investigated the transition 
from bound to scattering states. 

Furthermore we discussed a repulsive delta-shell potential as a simple model 
showing resonances. Resonances can still be identified in the nonlinear case, though 
the definition of a resonance becomes somewhat ambitious. Two major effects of 
the nonlinearity were analyzed in detail: Firstly, the resonance positions are shifted 
proportionally to the effective nonlinearity and the resonance width increases with 
ges- Secondly, scattering states cease to exist with an increasing repulsive nonlinearity, 
whereas resonances survive longest. 
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